Search for: All records

Abstract. Geologic archives of the Laurentide Ice Sheet (LIS) provide abundant constraints regarding the size and extent of the ice sheet during the Last Glacial Maximum (LGM) and throughout the deglaciation. Direct observations of LGM LIS thickness are non-existent, however, due to ice surface elevations likely exceeding those of even the tallest summits in the northeastern United States (NE USA). Geomorphic and isotopic data from mountains across the NE USA can inform basal conditions, including the presence of warm- or cold-based regimes, while covered by ice. While warm-based ice and erosive conditions likely existed on the flanks of these summits and throughout neighboring valleys, cosmogenic nuclide inheritance and frost-riven blockfields on summits suggest ineffective glacial erosion and cold-based ice conditions. Geologic reconstructions indicate that a complex erosional and thermal regime likely existed across the NE USA sometime during and after the LGM, although this has not been confirmed by ice sheet models. Instead, current ice sheet models simulate warm-based ice conditions across this region, with disagreement likely arising from the use of low-resolution meshes (e.g., > 20 km) which are unable to resolve the high bedrock relief across the NE USA that strongly influenced overall ice flow and the complex LIS thermal state. Here, we use a newer-generation ice sheet model, the Ice-sheet and Sea-level System Model (ISSM), to simulate the LGM conditions of the LIS across the NE USA and in three localities with high bedrock relief (Adirondack Mountains, White Mountains, and Mount Katahdin), with results confirming the existence of a complex thermal regime as interpreted from the geologic data. The model uses a small ensemble of LGM climate boundary conditions and a high-resolution horizontal mesh that resolves bedrock features down to 30 m to reconstruct LGM ice flow, ice thickness, and thermal conditions. These results indicate that, across the NE USA, polythermal conditions existed during the LGM. While the majority of this domain is simulated to be warm-based, cold-based ice persists where ice velocities are slow (< 15 m yr−1), particularly across regional ice divides (e.g., Adirondack Mountains). Additionally, sharp thermal boundaries are simulated where cold-based ice across high-elevation summits (White Mountains and Mount Katahdin) flanks warm-based ice in adjacent valleys. We find that the elevation of this simulated thermal boundary ranges between 800–1500 m, largely supporting geologic interpretations that polythermal ice conditions existed across the NE USA during the LGM; however, this boundary varies geographically. In general, we show that a model using a finer spatial resolution compared to older models can simulate the polythermal conditions captured in the geologic data, with the model output being of potential utility for site selection in future geologic studies and for geomorphic interpretations of landscape evolution. 

Abstract. The Greenland Ice Sheet's negative mass balance is driven by a sensitivity to both a warming atmosphere and ocean. The fidelity of ice-sheet models in accounting for ice-ocean interaction is inherently uncertain and often constrained against recent fluctuations in the ice-sheet margin from the previous decades. The geological record can be utilised to contextualise ice-sheet mass loss and understand the drivers of changes at the marine margin across climatic shifts and previous extended warm periods, aiding our understanding of future ice-sheet behaviour. Here, we use the Ice-sheet and Sea-level System Model (ISSM) to explore the Holocene evolution of Ryder Glacier draining into Sherard Osborn Fjord, Northern Greenland. Our modelling results are constrained with terrestrial reconstructions of the paleo-ice sheet margin and an extensive marine sediment record from Sherard Osborn Fjord that details ice dynamics over the past 12.5 ka years. By employing a consistent mesh resolution of <1 km at the ice-ocean boundary, we assess the importance of atmospheric and oceanic changes to Ryder Glacier's Holocene behaviour. Our simulations show that the initial retreat of the ice margin after the Younger Dryas cold period was driven by a warming climate and the resulting fluctuations in Surface Mass Balance. Changing atmospheric conditions remain the first order control in the timing of ice retreat during the Holocene. We find ice-ocean interactions become increasingly fundamental to Ryder's retreat in the mid-Holocene; with higher than contemporary melt rates required to force grounding line retreat and capture the collapse of the ice tongue during the Holocene Thermal Maximum. Regrowth of the tongue during the neo-glacial cooling of the late Holocene is necessary to advance both the terrestrial and marine margins of the glacier. Our results stress the importance of accurately resolving the ice-ocean interface in modelling efforts over centennial and millennial time scales, in particular the role of floating ice tongues and submarine melt, and provide vital analogous for the future evolution of Ryder in a warming climate. 

Abstract. Geologic evidence of the Laurentide Ice Sheet (LIS) provides abundant constraints on the areal extent of the ice sheet during the Last Glacial Maximum (LGM). Direct observations of LGM LIS thickness are non-existent, however, with most geologic data across high elevation summits in the Northeastern United States (NE USA) often showing signs of inheritance, indicative of weakly erosive ice flow and the presence of cold-based ice. While warm-based ice and erosive conditions likely existed on the flanks of these summits and throughout neighboring valleys, summit inheritance issues have hampered efforts to constrain the timing of the emergence of ice-free conditions at high elevation summits. These geomorphic reconstructions indicate that a complex erosional and thermal regime likely existed across the southeasternmost extent of the LIS sometime during the LGM, although this has not been confirmed by ice sheet models. Instead, current ice sheet models simulate warm-based ice conditions across this region, with disagreement likely arising from the use of low resolution meshes (e.g., >20 km) which are unable to resolve the high bedrock relief across this region that strongly influenced overall ice flow and the complex LIS thermal state. Here we use a newer generation ice sheet model, the Ice-sheet and Sea-level System Model (ISSM), to simulate the LGM conditions of the LIS across the NE USA and at 3 localities with high bedrock relief (Adirondack Mountains, White Mountains, and Mount Katahdin), with results confirming the existence of a complex thermal regime as interpreted by the geologic data. The model uses higher-order physics, a small ensemble of LGM climate boundary conditions, and a high-resolution horizontal mesh that resolves bedrock features down to 30 meters to reconstruct LGM ice flow, ice thickness, and thermal conditions. These results indicate that across the NE USA, polythermal conditions existed during the LGM. While the majority of this domain is simulated to be warm-based, cold-based ice persists where ice velocities are slow (<15 m/yr) particularly across regional ice divides (e.g., Adirondacks). Additionally, sharp thermal boundaries are simulated where cold-based ice across high elevation summits (White Mountains and Mount Katahdin) flank warm-based ice in adjacent valleys. Because geologic data is geographically limited, these high-resolution simulations can help fill gaps in our understanding of the geographical distribution of the polythermal ice during the LGM. We find that the elevation of this simulated thermal boundary ranges between 800–1500 meters, largely supporting geologic interpretations that polythermal ice conditions existed across NE USA during the LGM, however this boundary varies geographically. In general, we show that a model with finer spatial resolution and higher order physics is able to simulate the polythermal conditions captured in the geologic data, with model output being of potential utility for site selection in future geologic studies and geomorphic interpretation of landscape evolution. 

Abstract. Studying the retreat of the Patagonian Ice Sheet (PIS) during the last deglaciation represents an important opportunity to understand how ice sheets outside the polar regions have responded to deglacial changes in temperature and large-scale atmospheric circulation. At the northernmost extension of the PIS during the Last Glacial Maximum (LGM), the Chilean Lake District (CLD) was influenced by the southern westerly winds (SWW), which strongly modulated the hydrologic and heat budgets of the region. Despite progress in constraining the nature and timing of deglacial ice retreat across this area, considerable uncertainty in the glacial history still exists due to a lack of geologic constraints on past ice margin change. Where the glacial chronology is lacking, ice sheet models can provide important insight into our understanding of the characteristics and drivers of deglacial ice retreat. Here we apply the Ice Sheet and Sea-level System Model (ISSM) to simulate the LGM and last deglacial ice history of the PIS across the CLD at high spatial resolution (450 m). We present a transient simulation of ice margin change across the last deglaciation using climate inputs from the National Center for Atmospheric Research Community Climate System Model (CCSM3) Trace-21ka experiment. At the LGM, the simulated ice extent across the CLD agrees well with the most comprehensive reconstruction of PIS ice history (PATICE). Coincident with deglacial warming, ice retreat ensues after 19 ka, with large-scale ice retreat occurring across the CLD between 18 and 16.5 ka. By 17 ka, the northern portion of the CLD becomes ice free, and by 15 ka, ice only persists at high elevations as mountain glaciers and small ice caps. Our simulated ice history agrees well with PATICE for early deglacial ice retreat but diverges at and after 15 ka, where the geologic reconstruction suggests the persistence of an ice cap across the southern CLD until 10 ka. However, given the high uncertainty in the geologic reconstruction of the PIS across the CLD during the later deglaciation, this work emphasizes a need for improved geologic constraints on past ice margin change. While deglacial warming drove the ice retreat across this region, sensitivity tests reveal that modest variations in wintertime precipitation (∼10 %) can modulate the pacing of ice retreat by up to 2 ka, which has implications when comparing simulated outputs of ice margin change to geologic reconstructions. While we find that TraCE-21ka simulates large-scale changes in the SWW across the CLD that are consistent with regional paleoclimate reconstructions, the magnitude of the simulated precipitation changes is smaller than what is found in proxy records. From our sensitivity analysis, we can deduce that larger anomalies in precipitation, as found in paleoclimate proxies, may have had a large impact on modulating the magnitude and timing of deglacial ice retreat. This fact highlights an additional need for better constraints on the deglacial change in strength, position, and extent of the SWW as it relates to understanding the drivers of deglacial PIS behavior. 

Abstract. Numerical simulations of the Greenland Ice Sheet (GrIS) over geologictimescales can greatly improve our knowledge of the critical factors drivingGrIS demise during climatically warm periods, which has clear relevance forbetter predicting GrIS behavior over the upcoming centuries. To assess thefidelity of these modeling efforts, however, observational constraints ofpast ice sheet change are needed. Across southwestern Greenland, geologicrecords detail Holocene ice retreat across both terrestrial-based and marine-terminating environments, providing an ideal opportunity to rigorouslybenchmark model simulations against geologic reconstructions of ice sheetchange. Here, we present regional ice sheet modeling results using theIce-sheet and Sea-level System Model (ISSM) of Holocene ice sheet historyacross an extensive fjord region in southwestern Greenland covering thelandscape around the Kangiata Nunaata Sermia (KNS) glacier and extendingoutward along the 200 km Nuup Kangerula (Godthåbsfjord). Oursimulations, forced by reconstructions of Holocene climate and recentlyimplemented calving laws, assess the sensitivity of ice retreat across theKNS region to atmospheric and oceanic forcing. Our simulations reveal thatthe geologically reconstructed ice retreat across the terrestrial landscapein the study area was likely driven by fluctuations in surface mass balancein response to Early Holocene warming – and was likely not influencedsignificantly by the response of adjacent outlet glaciers to calving andocean-induced melting. The impact of ice calving within fjords, however,plays a significant role by enhancing ice discharge at the terminus, leadingto interior thinning up to the ice divide that is consistent withreconstructed magnitudes of Early Holocene ice thinning. Our results,benchmarked against geologic constraints of past ice-margin change, suggestthat while calving did not strongly influence Holocene ice-margin migrationacross terrestrial portions of the KNS forefield, it strongly impactedregional mass loss. While these results imply that the implementation andresolution of ice calving in paleo-ice-flow models is important towardsmaking more robust estimations of past ice mass change, they also illustratethe importance these processes have on contemporary and future long-term icemass change across similar fjord-dominated regions of the GrIS.